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Introduction 

Sri T^nka hosts a population of 21 million people in a relatively small area of 
65,610 km 2 (Fig. 1). Relatively high rainfall (1830 mm annually) and an equable 
temperature regime along with fertile soils have enabled the adequate provision of 
food and water to its people. Although agricultural production has increased over 
the last century, the population has been rising, reaching 21 million in 2011 from 
14.7 million in 1980 (Department of Census and Statistics, 2012b). While advances 
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m a^ituna research have been noteworthy, the challenge of securing food supply 

exa ^ r by conversion of agricultural lands, internal conflict, linSS 
technological inputs, and periodic shortages of water that constrain yield. 

in mSTST ““T 1 ? 1 U% rf 46 nati0Ml ^ ss domestic product (GDP) 
2012 while providing direct employment to 31% of the labor force (Central 

of the o "l . ^ 20B)> in&ectl y contributing to the livelihoods of up to 70% 
population. Among the crops that account for 78% of the agricultural GDP 

Z ^ ’ ^iZ C ° nUt Cr ° PS WWch COntribute l3A %> 10.1% and 10 9% to 

die agncuitural GDP, respectively. Rice is the principal food crop and is grown 
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water. Around 879,000 fanner families are engaged in rice cultivation (20% of the 

STd IZr 8 f'Sfrt ,OT Up *° 32 * of POP" 13 '™ (Department of 
Craaus and Statistics of Sn Lanka, 2012a). The country has become nearly self- 

sufficient in nee production (Department of Agriculture, 2009). 7 
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Sri T-anka has achieved tremendous progress since 1950 in crop production 
and food availability (Department of Agriculture, 2010; Ministry of Agriculture 
and Agrarian Services, 2008). Yields grew at an impressive rate until leveling off 
in the mid-eighties (Department of Agriculture, 2009). Sri Lanka’s population is 
anticipated to continue to grow in the coming decades (Department of Census 
and Statistics, 2012a), creating an ever-greater demand for food security on the 
household, sub-district, regional, and national scales (Sangakkara and Nissanka, 


2008). 

The agricultural sector in Sri Lanka is vulnerable to climate shocks (Yahiya 
et al., 2011; Zubair, 2002). An unusual succession of droughts and floods from 
2008 to 2014 has led to both booms and busts in agricultural production, which 
were reflected in food prices (Gadgilcf al., 2011; Zubaif, 2005). In both instances, 
the majority of farmers and consumers were adversely affected. The cultivation of 
some of the other crops (such as onions, chilies, and vegetables) has become more 
economically viable especially with the constraint of water. A typical small-scale 
fa rmin g system is represented schematically in Fig. 2. The household unit includes 
a home garden, other field crops, and livestock. 
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Commercial Fertilizer Field 

Fig. 2. A schematic of the rice-dominated agricultural system in this region 
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There is a millennium-long history of adaptation (spatially and seasonally) to 
climate variability across the island through choice of climatically and environmen- 
tally appropriate varieties, tailoring of planting season by region, development of 
appropriate irrigation infrastructure, and social and trade arrangements to suit com- 
munal agricultural practices (Zubair, 2005). These adaptation methods afford some 
potential for coping with future climate change, but they may not be capable of 
dealing with unprecedented changes. There is an ongoing program of rice varietal 
improvement at the Department of Agriculture and universities which historically 


have an excellent track record. 

The AgMIP (Rosenzweig etal, 2012) regional assessment in Sri Lanka addresses 
these issues by generating future climate projections and using them to drive crop 
models that feed into integrated climate-crop-economic models that in turn quantify 
climate impacts on livelihoods and poverty. We also test adaptation strategies in this 
modeling framework under future representative agricultural pathways (RAPs). 


Interactions with Stakeholders, Adaptation Options, and Agricultural 
Pathways 

Government policymakers, water-resource managers, local government officials, 
and farmers have highlighted the need for reliable climate change information, the 
need to orient research to address the questions they face for allocating resources, 
and the need for researchers to stay engaged and to disseminate learning. 
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concentrations rise. The 
the rice crop approaches 
2005). Rising sea levels 
> salination ( Minis try of 
>f the economic impacts 
nate projections to yield 
ulasuriya et at, 2007) 
n climate and yield pro- 
rise in temperature and 
tions from live GCMs. 
r has been limited with 
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Interactions with stakeholders 

Personnel from the Department of Agriculture involved with the AgMIP regional 
integrated assessment have solicited advice at departmental meetings and are in a 
position to contribute to departmental policy. 

We have reviewed climate risk information needs with personnel of the Mahaweli 
Authority of Sri Lanka (MASL), which is a large river basin management agency 
that provides services for the farmers in its domain (including water distribution, 
seed distribution, agricultural extension, community development, and enterprise 
development). The MASL coordinates national water management, to which effort 
the project investigators have provided a weekly climate advisory. 

The University of Ruhuna researchers are engaged in a project to undertake 
multi-sectoral, multi-institutional, community-based climate risk management in 
the Nilwala river basin in the south (Zubair et al., 2014). The perspectives of 
farmers and officials in this region regarding climate change and adaptation were 
solicited. 
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We also interacted with bulk purchasers of rice husks for energy conversion and 
rice for export as regards their needs for climate based predictions and projections 
of rice harvests (Zubair et al., 2014). 


Key decisions for national policymaking 

Among the key decisions for rice farming systems for the decades ahead are devel- 
opment of national plans for agriculture, land and water resources, food security, 
environment, poverty alleviation, and guidance for long-term infrastructural invest- 
ments. 



Potential adaptation strategies 

Based on a review of national policy documents related to agriculture and climate 
change adaptation (e.g., Ministiy of Agricultural Development and Agrarian Ser- 
vices, 2007; Ministry of Environment, 2010), assessments of researchers, and inter- 
actions with stakeholders including farmers, the following adaptation strategies have 
been proposed: (1) investment in research, innovation, and technology diffusion; (2) 
promotion of better-adapted cultivation practices and adaptive management of farm- 
ing systems; (3) risk management of shocks to the farm economy due to climate, 
pest, diseases, and market volatility; (4) decrease in food losses from farm gate to 
consumption; and (5) restructuring of policies on farming subsidies, price supports, 
crop insurance, and trading policies. These policy interventions can be sharpened 
and prioritized with better information on climate change and its impacts. 


Farmer perceptions of climate change and potential adaption 

Researchers solicited farmer perspectives in the Kurunegala (Migalewa, Kadawara- 
mulla), Ampara (Mabawanawella), and Matara districts (Nilwala Basin) (Rg. 1) — 
details in Zubair et al. (2014). Farmers in Migalewa perceive the following as major 
climate-related problems: 

• Climate: (1) Changes in total seasonal rainfall; (2) changes in starting date and 
duration (number of days) of growing season; (3) increase in occurrence of 
drought; (4) more frequent drought following intense rainfall. 

• Crops: Increase in occurrence of crop damages and reduction in harvest due to 
drought. 

• Water supply: (1) drop of groundwater level; (2) decrease of water supply for 
cultivation; (3) reduction in the duration of water supply for cultivation. 

• Pest and diseases: (1) appearance of new diseases and pest attacks. 
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While issues related to flooding and deforestation were not discounted, these 
were not a priority for the farmers. They identified the following as practical adap- 
tation options: (1) application of organic fertilizers, insecticides, and herbicides; (2) 
use of drought-tolerant varieties and shorter-duration crops; (3) change of seasons for 
cultivation and of planting dates; (4) crop diversification and intercropping; (5) rain- 
water harvesting; and (6) diversification of livelihoods and non-agricultural income 
sources. Farmers did not think it was practical either to undertake drip irrigation or 
energy-intensive cultivation. 

Constructing Representative Agricultural Pathways 

Representative agricultural pathways (RAPs) ensemble an overall narrative descrip- 
tion of a plausible future development pathway and kdy variables with qualitative 
storylines and quantitative trends — consistent with higher level pathways such 
as SSPs and RAPs developed by the AgMEP global modeling group (Rosenzweig, 
et al., 2013a ) described further in Part 1, Chapter 5 in this volume. 

RAPs were constructed in consultation with researchers and officials, officers 
of the Department of Agriculture, and the other departments based on literature 
from several government agencies, consultation meetings, and iterative review. The 
salient projected changes that are most likely in the coming decades are: 

9 The current trend of economic growth is projected to continue, 
o Population growth rate remains low and urban population increases while rural 
household size decreases, and farm size declines. 

. Rice productivity shall be adversely affected due to deterioration of soil fertility, 
extreme weather events, climate change, and phasing out of fertilizer subsidies 
(without adaptation). 

. Transportation infrastructure improves and trade policy encourages exports and 
imposes high taxes on imports of food that can be grown locally. 

• Technology developments in the rice sector with high-yielding, drought- and 
flood-tolerant cultivars, and improved crop management practices shall improve 
rice productivity. 

Summary Narrative of the Adopted RAP 

Best practices are implemented in a context of economic growth, low population 
growth, urbanization, and increasing productivity in the rice sector in a business- 
as-usual scenario. The government invests more in agriculture to improve food 
security through self-sufficiency in rice within a framework that promotes the ability 
of the rice sector to cope with impacts of a variable climate. This is done in a 
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magnitude of change and the confidence of all the variables were assessed with the 
involvement of experts. 

In relation to the biophysical parameters, temperature changes and land degrada- 
tion were identified as important, while there was low confidence on the precipitation 
increases that have been projected in the past. Among the socio-economic parame- 
ters, the current and likely future dearth of agricultural labor and outward migration 
for employment were identified as important parameters. Among the institutional 
and policy changes, a reduction in fertilizer subsidies is anticipated, even while 
investment in the agriculture sector could increase and the coordination among dif- 
ferent agencies dealing with agriculture could improve. In the technological realm, it 
is anticipated with high confidence that farmers will have improved access to infor- 
mation on markets, weather, and climate in the future and that improved varieties 
of rice will become available. 

Data 

Sri Lanka has rich archives of agronomic data (Pain, 1986; Wickremasinghe, 2006), 
which have led to (1) long-term studies on crops, including detailed field trials and 
coordinated rice trials; (2) long-term climate data at fine resolution (e.g., Zubair, 
2004); (3) household-level survey data and aggregated data for economic analysis 
and (4) data from rice research stations and university research programs. 


Climate data 

Station data 

Meteorological data for Batalagoda (7.51° N, 80°43' E), Maha Illuppallama (8° 10’ 
N, 80°45’ E), and Maho (7°82’ N, 80°27’ E) of precipitation, maximum temperature, 
and minimum temperature were obtained for the period of 1 980-2012. Daily rainfall, 
minimum and maximum temperature, atmospheric pressure, and sunshine hours are 
recorded at these stations (Zubair et al., 2014). 

Reanalyses data 

Hie AgMERRA data produced by the AgMIP Climate Team (Ruane et al., 2014) 
from the MERRA data (Rieneckcr et al, 2011) were used for gap-filling and for 
variables that were not available, after applying necessary bias-corrections. 

Climate change model simulations and scenarios 

The data from the 20 GCMs in the CMIP5 archive (Taylor et al., 2012) were 
used for projecting future climate. Analysis was restricted to five GCMs (CCSM4, 


/. Zubair et ad. 


GFDL-ESM2M, HadGEM2-ES, MIROC5, and MPI-ESM-MR) that were idenh- ^ 
lied by AgMP as relatively skillful in reproducing the climate over South Asia bv 
AgMIP. 1 


Crop Data '■£ 

Aggregate historical data on area of cultivation, production, and yield are estimated ■ ' 
collectively by several government agencies. "V 

• :* *? 

Rice experimental trials - 

Experimental data for rice trials are available from the Department of Agriculture ? 
(Batalagoda, Maha Illupallama, and other field stations). Calibrations were carried £ 
out on the following varieties of data in wide spread use in the region — Bg300 4 
(three months duration), At308 (three months), Bg357 (three-and-a-half months), £ 
Bg358 (three-and-a-half months), and Bg379-2 (four to four-and-a-half months)! ' 
Further details are available in the chapters on calibration in Zubair et al. (2014) 
including for other varieties (Bg250, Bg304, Bg307, Bg94-1). 

Data from (1) detailed field trails and (2) Coordinated Rice Variety Trials & 
(CRVTs), were used for caHbration and evaluation of models. These detailed exper- ft 
iments, conducted at the Rice Research and Development Institute (RRDI) of Sn 
Lanka at Bathalagoda, from 2001-2010 in the major (Maha) and minor (Yala) 
growing seasons with two nitrogen fertilizer levels. CRVTs were conducted at two 
different agroecological regions of RRDI and FCRDI (Field Crop Research and 
Development Institute (FCRDI) at Maha llluppallama), during the two growing 
seasons from 2000 to 2008. 

The size of the experimental plot was 1 8 m 2 (6 m x 3 m). Recommended standard 
procedures by the Department of Agriculture (DOA) were adapted for rice crop 
management. Irrigation was provided, then the land was ploughed, crop residues of 
the previous season were incorporated three weeks before planting, and the soil was 
puddled continuously with subsequent ploughing/harrowing actions undertaken at 
least twice to make the fields weed-free and level before planting. 

The crop was established by broadcasting pre-germinated seeds (direct seeding 
at a rate of around 350 plants/m 2 ). Irrigation water was provided as required and 
maintained at a standing water level of about 25 mm starting from one week after 
emergence until grain-filling was completed. 

Fertilizer application was carried out as recommended by the DOA. Split applica- 
tion was practiced as atplanting, first top dressing (two weeks after planting), second 
top dressing (five weeks after planting), third top dressing (seven weeks after plant- 
ing) for three-month varieties, and additional application of N at booting stage for 
three-and-a-half and four month varieties. Medium fertilizer treatment represents 
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the DOA recommended fertilizer levels. Continuous monitoring of pest and disease 
outbreaks was undertaken and appropriate control measures were applied whenever 

^orthe detailed experiments, growth and development parameters (plant height, 
leaf area, and tiller count) were recorded at regular intervals, along with time taken 
to panicle initiation, anthesis and maturity stages. Yield and yield components were 
measured at maturity. For CRVT experiments, time taken to heading and harvesting 
stages, and yield measurements were recorded. 

i>ata for Integrated Assessment 

Researcher-conducted surveys of fanning systems are available for a cross-section 
of farmers in selected locales in the Kurunegala district (Fig- 
referred to as Set 1 (Nikaweratiya, Rajangane, and Batalagoda) and Set 2 (Kada ara 
mulla, Migalewa) and details are provided in Zubair etal. (2014). All of the lotions 
except Kadawaramulla practice cultivation in the minor' Yah a s^omto Ka^wam- 
mulla there is a system of irrigation in alternate years for Yala for each bank of toe 
feeder irrigation canal. Including Kadawaramulla ensures that we capture practices 
that are representative of the entire district. Although we have conducted surveys 
for two successive Yala seasons for Kadawaramulla, the year it received negation 
coincided with a drought that devastated the rice crop during the Yala 
Set 2, to be able to represent both seasons, we reporthere on the results for Migalewa 
while detailed results for Kadawaramulla are not shown here, but axe included in 

Zubair et al. (2014) and in Herath et al. (2013). , . 

Soil and climate files were developed for the respective sites based on die 
recorded information (Mapa er al, 2005). Additional soil testing was undertaken for 
Kadawaramulla and Migalewa (Zubair er al, 2014). Crop management practices for 
farmer fields were estimated based on the farmer-survey information and standard 

practices. 


Methods of Study 


Climate 

The work followed methods outlined in the AgMIP Handbook (Rosenzweig * al., 
2013a) and described further in Part 1, Chapter 3 in this volume. 


Quality control 

Inconsistent data entries were identified and removed. Neighb oring station data were 
compared in cases of unusual observations. The solar radiation data were taken from 
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AgMERRA datasets developed by the AgMIP climate team (Ruane eta!., 2014) and 
compared with estimations from sunshine hours based on a previously calibrated 
algorithm (Samuel, 1991). 

Missing values of variables and gapsfrom data removed for quality control we«3 
fiUed in. If the gap was less than two days, it was filled based on the values on l 
both sides of the gap and after cross-checking with observed values of neighboring * 
stations for validation in the case of temperature. Other gaps of data (more than two ' 
days) were filled by using the AgMERRA dataset prepared by the AgMIP climate ✓ 
team with bias-correction. 

. In rasta ” ces where toe farm locations were not close to the meteorological sta- • 
tions, a method of extrapolation developed by the AgMIP Climate Team based on the : - 
high-resolution monthly spatial climatologies developed by WorldClim was used, f 

Mid-20th century projections :'i- 

^ « GCM based to to mid-20th centmj | 

UU4O-2U70) under a high-emission-concentration pathway (RCP8.5) (Moss et a! 

2010; Taylor et al., 2012). The results from the following five GCMs were used for I 

^ e “° re detaiJed analysis 35 re P resent toe South Asian region better- V 
C CSM4,GFDL-ESM2M, HadGEM2-ES, MIROC5, and MPI-ESM-MR (labeled f 
as E, I, K, O, and R respectively). The methodology for downscaling the climate ?■ 
projections to each station was based on the “delta method” (Wilby et al 2004) The * 
implementation is described in detail in Part 1 , Chapter 3 in this volume. Time senes 4 
at individual stations were created by modifying the meteorological station data with ¥ 
memi monthly spatial differences provided by the WorldClim dataset (Hijmans eta}., 4 

nf ° f " 1,118 3Verage ^^toro precipitation and 

of afi 20 GCMs were tested based on whether there was a change of two standard 
deviations of yearly variation (see Part 1, Chapter 3 in this volume). 4 

V 

Crop Model Calibration and Simulations 

Calibration "jf- 

Cdjhati°n of DSSAT-CERES-Rice version 4.5 (Jones et al. 2003), and APSIM- H 
ORYZA version 7.5 (Keating « al. 2003) crop models was undertaken for the 
commonly used varieties of rice in Sri Lanka for which experimental data wete V 
available (Bg. 50. Bg 300, Bg 357, Bg 358, Bg 379-2, At 307, and At 308). These if: 
data are based at Batalagoda and Maha Illupallama, along with the data from the f 
national CRVT for the same locations. ~ 


>. - 
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Table 2. Genetic coefficients generated for DSSAT-CERES rice crop model for four cultivars. 


PI P2R P5 


Bg300 

407.9 

134.7 

396.2 

Bg357 

382 

122.5 

393.8 

Bg 358 

435.6 

72.5 

500.1 

At 308 

400.8 

50.3 

400.7 

Bg 379-2 

662.9 

140.9 

437.1 


P20 

Gl 

G2 

G3 

G4 

12.6 

50.7 

0.026 

0.95 

1.1 

11.5 

55.5 

0.025 

0.95 

1.23 

10.7 

50.3 

0.019 

0.69 

1.1 

12 

68.5 

0.019 

1 

1.2 

10.8 

58.9 

0.03 

0.5 

1.16 


DSSAT CERES-Rice calibration 

Detailed experiments and CRVT from two experimental sites at Batalagoda and 
Mahailluppallama were used for generalized likelihood uncertainty estimation 
(GLUE) runs, initially only for phenology. Then GLUE-estimated phenology was 
fixed for the estimation of growth parameters (Gl, G2, G3, G4) by using measured 
growth variables in detailed experiments. The grain size (G2) was calculated for 
selected varieties based on historical data and fixed for the variety. Therefore, only 
Gl, G3, and G4 were estimated from GLUE runs of detailed experiments. The crop 
model simulations were evaluated by using model efficiency criteria for leaf-area 
index (LAI), tops weight, tiller number, and yield predictions for detailed experi- 
ments, and yield only for CRVTs. The genetic coefficients obtained are provided in 
Table 2. 


APSIM calibration 


The same dataset used for DSSAT calibration was used for APSIM calibration. The 
APSEM-ORYZA model parameters of rice according to different growing stages 
were considered for the calibration of the respective varieties; these include two 
photoperiod-sensitive stages and two growth parameters. Developmental stages of 
the crop model are: 


1 . DVLR; Development stage in juvenile phase 

2. DVRI; Development stage in photoperiod-sensitive j 
phase 

3 . DVRP; Development rate in panicle 1 
development 

4. DVRR; Development rate in 
reproductive phase 

The fixed values of DVS at important growth stages in APSIM are 0 (at sowing), 


Photoperiod-sensitive 


Photoperiod-insensitive parameters 


0.65 (at panicle initiation), 1 (at flowering), and 2 (at maturity). 

Developmental rate parameters (DVS) were calculated by using growing degree 
days (GDD). GDD were calculated by using T max , T min , and base temperature (T b ), 
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Table 3. Genetic coefficients generated for APSIM-Oryza rice crop model for four cultivars 


Model 

Bg 300 

Bg 357 

Bg 358 

Bg 379-2 

parameter 

Yala 

Maha 

Yala 

Maha 

Yala 

Maha 


DVRJ 

DVRI 

DVRP 

DVRR 

MOPP 

0.001071 

0.000853 

0.000646 

0.001755 

13 

0.001071 

0.001071 

0.000646 

0.001755 

13 

0.000878 

0.000737 

0.000636 

0.001729 

13 

0.000878 

0.000878 

0.000636 

0.001729 

13 

0.000869 

0.000721 

0.000641 

0.001871 

13 

0.000869 
0.000869 
0.000641 
0.001871 . 
13 

0.001037 V 
0.000285 * 
0.000643 V 
0.001869 v 
11 5 


i.e., 8 C. GDD for phenologica] events were calculated based on the ORYZA-2000 
model. The parameters estimated for APSIM inputs are given in Table 3. 




DSSAT-CERES and APSIM-ORYZA crop model simulations 

DSSAT-CERES and APSIM-ORYZA simulations were undertaken based on the 
former surveys identified as Set 1 and Set 2. These simulations were undertaken with 
consistent genetic coefficients, consistent choices for management but differentiated 
by the details of the farm surveys and the soil and weather for Set 1 and Set 2. 

The “ rr tools ” Sloped by AgMIP were used to undertake the large number of - i 
simulations required in a batch mode. The weather, soil, fertilizer, water, and crop S \ 
management data were incorporated into the input files (i.e., field overlay file, survey 
data file, and seasonal strategy file) that were prepared for the QUADUI AgMIP IT •; 

tool (Rosenzweig etal., 2013a). This enabled the generation of DOME and ACMO 
files for input into the crop simulations and integrated assessments. ''•/ 

Farm-survey data were used in order to match crop management (Hoogenboom ' : 
et al., 2010), i.e., planting, harvesting, application of inorganic fertilizer, irrigation, f ; 
and applications of crop residues and organic manure. 

Yields for the “current year” (2012/2013 for Set 1), were simulated for individual . • 
farm fields. For Set 2, the current year was September 2011 to August 2012 for 
Kadawaramulla and September 2012 to August 2013 for Migalewa. 'v 

Using these simulated results for the “matched data case” (Rosenzweig et aL ^ 
2013 a) were obtained where the models use weather, soil, and management for each 
field to simulate production matched with observed yields for each field. Yields for 
the future periods (2040-2069) were obtained by using the corresponding climate 
change scenarios derived from global climate models and observed weather data 
(Rosenzweig etal., 2013a, 2013b). 

hi simulating the future period, automatic planting was set within an interval 
of time with conditions for soil water content (i.e., target rainfall accumulation) 
specified at 20 mm. The harvesting option was set as “when the crop is mature”. 
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In simulating the future period for Migalewa and Kadawaramulla, automatic 
planting was set for one-day duration for the planting window and for soil water 
content specified at 150 mm. 


Integrated Assessment 

We used the Tradeoff Analysis Model for Multi-dimensional Impact Assessment 
(TOA-MD) economic model, which is a multi-dimensional impact assessment pro- 
gram that utilizes statistics of farm characteristics (means, standard deviations, coef- 
ficients of variation) to simulate the impacts of environmental changes on economic 
conditions (Antle, 2011; Antle and Valdivia, 2006). In this study we use it to assess 
climate change impacts on socio-economic aspects of farming systems with the 
intent of providing, support for informed policymaking. The model could be used in 
the future for prospective evaluation of technology adoption and technology impact 
assessment, payments for ecosystem services and environmental change impact and 
adaptation. 

The TOA-MD analysis was undertaken for a population of heterogenous farms 
under one stratum with one production activity of growing rice. These methods were 
followed to answer the following three “core climate impact questions” as posed by 
AgMIP (Rosehzweig et al, 2013a). 

1. What is the sensitivity of current agricultural production systems to climate 
change? 

This seeks to address the impact of climate change on the current production. 

2. What is the impact of climate change on future agricultural production systems? 
This seeks to address the impact of climate change on the future production 
system. Technological and other developments in the agricultural sector may 
alter the current production system in the future. Hence vulnerability to climate 
change may be dissimilar to the case of the current production system, which is 
described in the Question 1 simulations. 

3. What are the benefits of climate change adaptations? 

This line of inquiry seeks to assess the advantages of adaptation strategies in 
abating deleterious climate change impacts on the future rice-production system. 
While we have undertaken work on adaptation, it shall be reported in the future. 

While TOA-MD analysis was initiated for all of the farming systems, in this 
chapter results are presented only for Migalewa. 

In setting up the TOA-MD model there was one key assumption. The farmer 
surveys in Migalewa were undertaken separately for each of the two seasons in 
the same small region. The populations of farmers for each season have similar 
household economics and use similar crop management practices. In fact, nine 
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Table 4. 


Farm sizes and household sizes for the two surveys taken durine 
Vlaha seasons conducted at Mi palewa g 


Season 

Sample size 

Farm size 

CV of farm size 

HHsize 

CV of HH 

Maha 

Yala 

31 

31 

0.89 

0.82 

22.16 

33.95 

4.16 

4.45 

29.17 

28.91 


- 


System H of the Mahaweli River Basin AaK£ £ ** 

should have similar practices, market conditions, technical resources and 
histones (Raby and Merrey, 1989). Based on this assumption the results fr™ S 
seasons were combined to provide annual estimates. ° m ^ 

The two seasons were considered as two crop activities in the TOA un M ^ , 
setup. Two sets of crop simulations were undertaken to address the mrP - 

w,«h 0* use of ciireaie project from toe five le^Ms “ 

these two crop simulations were for: ’ Y 

1- Current climate with current production system Y 

2. Future climate with current production system 

^s=5aS3S5S 

specified for Sri Lanka (Grom 1 l ^ ^ " h » 

Results 
Climate results 

The climate projections for the mid-century (2040-2069^ for n„t Q i a 
on 20 GCMs uader higher greenhouse gas !o.e^„ 
ftg. 3. There is a coasisien, rise in temperature in “ 

The rucre.se iu temperature is near* consist across 2 «or “a 




V- 


* 
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Tlble 5. Values for trends for variable, and Moors ™ed in the TOA-MD that were obnnnnd 
mostly from Groups 1 and 2. — _ — — 


Group 1-Trend 
factors mostly 
from RAPs 


Group 2-Trend 
factors mostly from 
IMPACT global model 

Vaiable 

Description 

Value 

Description Value 

Yield trend x 

RAP defined value 

1.25 

AgMIP Reference scenario 1 498 

under IMPACT Global 


Price trend <p 1 


Price trend tp2 


Household size 
(HH> 

Fann area (Ha) 
Cost trend (40 


Non- 

agricultural 
income (Rs) 


AgMIP Reference price 
under IMPACT model 

RAP defined value 
(assumed to be matched 
with cost trend) 

RAP defined value 
(percentage decrease) 
RAP defined value 
(percentage decrease) 
RAP defined value 


RAP defined value 
(percentage increase) 


model 

1.18 -AgMIP Reference price trend 
under IMPACT Global 
model ■ 

1.2 SSP2_HGEM-DSSAT_5crop 
scenario price Trend under 
Impact Global model 
- 16% RAP defined value 

(percentage decrease) 

— 18% RAP defined value 

(percentage decrease) 

1 .2 Cost from IMPACT model 
(assumed to be same as 
SSP2 JHGEMLDSS AT _5crop 
scenario price trend ) 

30% RAP defined value 

(percentage increase) 


1.179 

1.462 

-16% 

-18% 

1.462 

30% 


season, there is an increase in temperature ranging from 1.5 to 2.8°C »d for fire 
Maha season there is an increase in temperature ranging from 1.1 to 2.4 C. Th 
monthly climatology of rainfall in the future period shows a significant mciease 
ftom September to November for all GCM scenarios and less change m other 

m ° f significance test was undertaken with a change of two standard deviations of 
yeariy Nation for either parameter taken as significant for the mean alone, and 
S mlT and variability-based scenarios (see Part 1, Chapter 3 m thrs volume k 
The* results show that the temperature rise is significant in both seenanos for al 

penrid (with one GCM showing a drop of 0.1 mm/day) from a current avemgerf 6 
L/day. This increase in rainfall is significant in bo* tearing scenarios fo befit Yala 
Zl Maha for 16 of the 20 GCM outputs including four of the five selected GC 
(except for HADGEM2-ES). Note this exceptional GCM was one of the warm . 
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Fig. 3. Top left: The monthly and seasonal climatology for temperature in the present (1980-2010 
as a solid line) and m the future projections (2040-2070 as box-and-whisker plots for the 20 GCMs) 
or atalagoda in northwest Sn Lanka. The line shows the climatology for one GCM (COM4) as an 

GCMS md th ° Se Sdccted for detailed analysis are shown 
in red. These are the GCMs CCSM4. GFDL-ESM2M, HadGEM2-ES, MIROC5, and MPI-ESM-MR 

an are labeled as E, I, K, O, and R. Top right: Shows the same but at the annual time-step (“ann”) and 
or toe scions of Januaty to March (JEM), April to June (AMI), July to September (JAS) and Octobe. 

. , u « em ber(ONp>. rhese seasons approximately correspond to early and late Maha (Octobei to 
March) and Yala (Apnl to September) agricultural seasons. Second panel, left and right: The same as 
at the top but for precipitation rather than temperature. Third panel: The bottom graphic shows the 
seasonally averaged current (solid square) and future projections (letters as in top panel) and climate 
projections or the mid-ccntuiy (2040-2069) under RCP8.5 for Batalagoda. The months in the season 
most hkcly to affect the nee crop are chosen - October to February for Maha and April to July to. 


1H 
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agoda, Sri Lanka) 



j?M Aifij ~MS Olto 



For the Yala season, there is an increase in rainfall of up to 2 mm/day for 10 of 
the GCMs and a decrease in rainfall for the rest of up to 1 mm/day from an average 
of 3.5 mm/ day However, only 10% of the GCMs showed a significant anomaly (a 
rise in this instance), and none of the five selected GCMs had a significant change 
in rainfall. 


goda, Sri Lanka) 



aim jHT auj ^Xs~ om~ 

muy RCP8 5 GCM. (Balalagoda, Sri Lanka) 


°« j* 


1 30 31 3 9. S3 

(M Ttmpermm fey 

e in the present (1980-2010 
sker plots for the 20 GCMs) 
"or one GCM (COM4) as an 
detailed analysis are shown 
IROC5, andMPI-ESM-MR 
innual time-step (“ann”) and 
iptember (JAS) and October 
and late Maha (October to 
left and right: The same as 
• bottom graphic shows the 
as in top panel) and climate 
a. The months in the season 
Maha and April to July for 


Crop modeling results 
DSSAT-CERES calibration results 

Simulated results with DSSAT-CERES for time taken to panicle initiation, anthesis, 
and maturity from emergence compared well with observed values for detailed and 
CRVTs experiments. RMSE (root mean square error) was around four to five days 
for all varieties (see Fig. 4 for calibration results of anthesis and yield for the Bg 300 
cultivar). RMSE for yield was lower for the Bg 357 cultivar (425 kg/ha), compared 
to the other two varieties Bg 300 and Bg 358, which ranged from 1800 to 2000 
kg/ha. All varieties had a similar RMSE of ~1400 for grain yield. Deviations for 
top weights and grain yields for these varieties compared to those observed were 
mainly due to biotic stress. For example, CRVTs of Bg 358 recorded lower observed 
yields due to a pest attack at grain-filling, which could not be controlled in time. 


APSIM-ORYZA calibration results 

The APSIM-ORYZA model calibration results (Fig. 4) showed a good fit, indicating 
RMSE for time taken to panicle initiation, anthesis, and maturity, were in the range 
of three to six days. RMSE for grain yield was around 1000 kg for Bg 300 and Bg 
357, while it was higher (1200 kg) for Bg 358. Deviation of simulation for grain 
yields for these varieties compared to the observed was mainly associated with the 
biotic stress. For example, CRVTs of Bg 358 recorded lower observed yields due to 
a pest attack at grain-filling, which could not be controlled in time. 


Crop simulation results for farm survey yields and climate scenario effects 

Parallel studies were undertaken for two sets of locales chosen from the Kurunegala 
district as described in the methodology. Set 1 included farm survey data from Bata- 
lagoda, Nikaweratiya, and Rajanganaya and Set 2 included Migalewa and Kadawara- 
mulla. 

Crop simulations were generated using DSSAT and APSIM crop models. The 
same genetic coefficients as described above were used. Simulations for the base 
year (2012-2013), historical period (1980-2010), and five GCMs at mid-century 




06s«v«l<bys CtsavidytelflWh,, 

Fig. 4. Comparison between simulated and observed yields with a 1:1 line (— ) of (A) days to ^ 
anthesis and (B) grain yield for the rice variety Bg 300 using DSSAT (top row) and APSIM (bottom . 


(2040-2070) for the major and minor seasons were obtained. These results were •' 
calibrated and validated accordingly. :i - • 

SET 1 : Batalagoda, Nikaweratiya, and Rajanganaya ';y 

Simulation results for the first set (collection of all three locations) for the major k' 
and minor seasons of the base year (2012/2013) are presented in Fig. 5. Signifi- *■' 
cant positive relationships were found between observed and simulated rice yields ‘ 
for both seasons for the DSSAT and APSIM crop models. A seasonal comparison £ 
revealed slightly greater overestimation for the minor season by DSSAT (RMSE of c - ■; 
1300 kg/ha) than for the major season (RMSE of 1200 kg/ha). 
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DSSAT-Major season: Kurunegala 2012/13 BSSAT-Mhor season: Kurunegala 2012/13 
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•cations) for the major 
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APSM-Major season. Kurunegala 2012/13 
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APSM-Minor season: Kurunegala 2012/13 
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/ 
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Pig 5 Comparison between DSSAT-CERES simulated and observed yields of farmer fields in 
Kurunegala (Set 1) for the base year of (20 12-201 3) with a 1 : 1 line (- - -). Thkis shown for simulations 
from DSSAT in the top panel for the major (Al) and minor (Bl) s^ons/Ihe reside from APSIM- 
ORYZA simulations are shown in the bottom panel for the major (A2) and minor (B2) seasons. 

However, APSIM simulated higher yields during the major season (RMSE of 
1700 kg/ha) compared to the minor season (Fig. 5). Multi-model comparison of the 
probability of exceedance of simulated vs. observed yield revealed a more-or-Jess 
similar relationship between the two crop models for the minor season compared to 
the major season (Fig. 5). 

Our simulations are undertaken with the assumption that the crop is not water- 
stressed. Water limitation for the rice crop is relatively lower in the major season 
than in the minor season. Sensitivity of these crop models to water status may have 
contributed to these differences. Farmer survey revealed that some fanners were 
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(a) Major season: Kurunegala 2012/13 



Observed yield fkfl/ha) 


(b) Minor season: Kurunegala 2012/13 



Observed yield (k<j/ha) : ■* 

Fig. 6. The probability of exceedance yields for observed and simulated (DSSAT and APS1M) yield = 
for farmer fields in Kurunegala for the base year (2012-2013) for the major (A — top) and the minoi 4 . 
(B — bottom) seasons. 

unable to get the required fertilizer at the correct time, and the application was ' 
delayed particularly in the minor season. Few farmers, especially in the Batalagoda - > 
region, also reported that their fields were infested with weeds that they were unable . I 
to control. These reasons may have also contributed to lower observed yield for 1 
some farmers, though the amount of fertilizer applied remained the same, resulting $ 
in higher simulated yields, RMSE, and poorer yield distribution fit for probability 
of exceedance in both crop models (Fig. 6). 

Yield prediction with DSSAT with GCM scenarios 

DSSAT results revealed that rice yield is reduced for all GCMs with a different mag- 
nitude for each season (Fig. 7). Compared to the baseline period, percentage yield 
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(a) DSSAT yield variations for 
major season in Korunegala 



(a) APSIM yield variations for 
major season in Kurunegala 



(b) DSSAT yield variations for 
minor season in Kurunegala 



M , , , - 

. BASEUNE E I K O S 

(b) APSIM yield variations 
for minor season in Kurunegala 



BASEL INC E I K 0 R 


Fig- 7. Top: DSSAT simulated rice yield for the major (A) and minor (B) seasons for five GCMs 
(mid-century) compared to baseline period for farmer fields. Middle: As in top panel for APSIM 
simulated rice yield for the major (A) and minor (B) seasons. Bottom: As top for DSSAT simulated 
rice yield for the major (A) and minor (B). 


d (DSSAT and APSIM) yield 
jor (A — top) and the minor 


id the application was 
.dally in the Batalagoda 
ts that they were unable 
wer observed yield for 
aed the same, resulting 
ution fit for probability 


Is with a different mag- 
;riod, percentage yield 


reductions for the five GCMs CCSM4, GFDL-ESM2M, HadGEM2-ES, MER.OC5, 
and MP1-ESM-MR (labelled E, 1, K, O, and R) were 13, 12, 19, 15, and 16%, while 
percent reduction for the minor season was 27, 28, 41, 33, and 36%. The highest 
reduction was observed for the HadGEM2-ES (K) GCM for the minor season being 
the warmest of selected five GCMs. Since the Yala is the warmer season, it is likely 
that if temperature rises excessively, then the rice crop which is usually at its higher 
range of temperature tolerance (Weerakoon et al, 2005) shall be adversely affected. 

Yield prediction with APSIM with GCM scenarios 

Different responses were observed for the climate change yield simulations with 
APSIM. Though a slight yield reduction was observed for all five GCMs in the 
major season, higher yields were simulated for four GCMs (except for the Had- 
GEM GCM) in the minor season. Percentage yield reductions observed for the 
Maha season compared to the base period for the five GCMs (CCSM4, GFDL- 
ESM2M, HadGEM2-ES, MIROC5, and MPI-ESM-MR) were 5, 5, 8, 5, and 6% 
respectively while only 2% reduction was observed for the HadGEM2-ES GCM for 
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the Yala season. The other GCMs (E, l O, and R) simulated 8%, 9%, 9%, and 3% i 
greater yields respectively for the minor season, compared to the baseline period. ' 
In both seasons the highest y ield reduction was observed for the HadGEM GCM. ‘r 

SET 2: Kadawaramulla and Migalewa v 

Farmer-specific rice yields were simulated for Set 2 for the historical period and I. 
the projections of the five GCMs of mid-century (2050s) for the major and minor f 
seasons using DSSAT and APSIM. :/ 

In the case of Kadawaramulla of Set 2, the rice cultivation was only available for 
the major season due to a devastating drought in the minor Yala season for which v : " 
the field survey results are available. We were not able to obtain data for the minor '% 
season in years before or after as there is a system of irrigation being provided only f ■ 
on alternate years for this region. The results for Kadawarmulla are included in % 
Zubair et al. (20 1 4). $ 

Results of the DSSAT simulations for the historical period of the major and : 'V' 
minor seasons for the Migalewa location of Set 2 (Figs. 8 and 9) are qualitatively 
similar to that presented for Set 1. The RMSE of simulated yields is mostly the v 
same (Fig. 8). The yield distribution as probability of exceedance matches well to 
the observed farm yield distribution (Fig. 9). T ■ 

The yields simulated with DSSAT are reduced for all five GCMs for both major % 
and minor seasons (Fig. 10). For Migalewa, the percentage yield losses in the H 
major season compared to baseline for the five GCMs (CCSM4, GFDL-ESM 2 M, A 
HadGEM2-ES, MIROC5, and MPI-ESM-MR) were 6.8, 6.1, 14.5, 7.5, and 7.9% ir- 
respectively. Similarly for the minor season, the yield reduction was 24.8, 26.7, 36.8 T 
24.6, and 30.6% respectively. For Kadawaramulla yield reductions in the major 
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Fig. 9. The probability of exceedance of observed and simulated yields with DSS AT during the base 
year (2012-2013) in Migalewa for major and minor seasons. 


(a) DSSAT Rice Yield Variation for 
Major Season in Migalewa 


(b) DSSAT Rice Yield Variation for 
Minor Season in Migalewa 


l 

£ 




Fig. 1 0. DSSAT simulated rice yield for the major (A) and minor (B) seasons for five GCMs (mid- 
centnry) compared to baseline period for farmer fields Migalewa. 

season were 3.1, 2.7, 9.8, 4.0, and 4.8%. The highest yield reductions were for 
simulations with GCM-HadGEM2-ES (Fig. 10). 


Integrated Assessment Results 

The results under Migalewa of Set 2 are presented under the integrated assessment 
section. The results from the TOA-MD analysis for Kadawaramulla are reported in 
Zubair et al. (2014) and Herath et al. (2013). 


Addressing Core Question l 

With current agricultural systems and future climate: The projected climate change 
leads to a negative economic impact on 75—85% of fanners in Migalewa (Table 6). 
The percent of gains under all GCMs show similar results of approximately 13% of 
mean net returns. Net returns per farm decrease by 14-26% depending on the GCM 
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Table 6. TOA-MD analysis showing yields, net returns, per capita income, and poverty rate *i«i 
and without climate change (top half labeled Question 1) and with climate change and changed 
agricultural practices (bottom half labeled Question 2) for Migalewa. The results are based on cl mL ' - 

change as assessed from five selected GCMsCCSM4,GFDL-ESM2M,HadGEM2-ES,MIROC5 and ’ 

MPI-ESM-MR (labeled B, I, K, O, and R) and trends as in Group 1, .Ji ' 



Stratum 

E 

I 

GCM 

K 

O 

• V 

- 

R ~S' ; 


Observed mean yield (Rice-Activity l-Maha 
season) (Kg/Ha/Season) 

4985 

4985 

4985 

4985 

4985 -y" 


Mean yield change (Rice) (%) 

-7 

-6 

-15 

-8 

9 p' 


Observed mean yield (Rice-Activity2-Ki/a 
season) (Kg/Ha/Season) 

■ 3508 

3508 

3508 

3508 

3508 


Mean yield change (Rice) (%) 

-25 

-b 

-37 

-25 

1 / v. - 

-31 ' 


Losers (%) 

76 

78 

85 

77 

80 '-v 


Gains (% mean net returns) 

13 

12 

13 

12 

13 

§ 

Losses (% mean net returns) 

-24 

-25 

-32 

-24 

-27 

l 

Observed net returns without climate change 
(Rs/Farm/Year) 

115777 

115777 

115777 

115777 

115777 


Observed net returns with climate change 
(Rs/Farm/Year) 

98630 

96632 

86290 

97697 

94200 ;S. 


Observed per capita income without climate 
change (Rs/Person/Year) 

54256 

54256 

54256 

54256 

54256 


Observed per capita income with climate 
change (Rs/Person/Year) 

49939 

49436 

46833 

49704 

48824 Vy.. 


Observed poverty rate without climate change 
(%) 

55 

55 

55 

55 

55 •l.» Vi 
' ■+ . 


Observed poverty rate with climate change 
{%) 

60 

60 

63 

60 

6i m y$ 


Projected mean yield (Rice- Activity l-Maha 
season) (Kg/Ha/Season) 

5830 

5830 

5830 

5830 

■ 

5830 


Mean yield change (Rice) (%) 

-7 

-6 

-15 

-8 

-8 •■r-'- 


Projected mean yield (Rice-Activity2- Yala 
season) (Kg/Ha/Season) 

5399 

5399 

5399 

5399 

5399 i f| 


Mean yield change (Rice) (%) 

-25 

-27 

-37 

-25 

-31 


Losers (%) 

77 

79 

88 

75 

82 >• 

CM 

Gains (% mean net returns) 

10 

10 

10 
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projection. The HADGCM-ES2 (identified by K) leads to the most significant losses 
and net impacts. Poverty rates increases by 8-15% for the five GCMs. Poverty line 
was defined as those earning less than $1.25/day. Highest poverty rate increase is 
shown under HADGCM-ES2 which has the highest temperature increase and the 
most yield drop. Note that the same rice prices are used for both System 1 and 2 with 
no trends used for yields and cost of cultivation. Hence analysis under Question 1 
gives a basic assessment of the socio-economic impacts due to climate change. 
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Addressing Core Question 2 

The impact of climate change on the future production system was assessed assuming 
that future production system will differ from current production system due to 
trends as described by the RAPs and the IMPACT global model and summarized 
in Group 1 in Table 5. The RAPs include an expectation of continuation of current 
development trends in the agriculture sector. The details of gains, losses, and net 
impacts are in Table 6 and in Fig. 1 1 . We also assessed the sensitivity to the choices 
of RAPs by using trend factors (Group 2) taken largely from the IMPACT model. 
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Fig. 1 1 . The gains, losses and net impacts as a percentage of mean farm net returns for Migalewa 
as estimated by the TOA-MD analysis for Question 1 ( top panel) of anticipated mid-century climate 
change and the current production system and Question 2 of anticipated mid-century climate change 
and future agricultural production system. The results are based on climate change as assessed inputs 
from five selected GCMs CCSM4, GFDL-ESM2M, HadGEM2-ES, MIROC5, and MPJ-ESM-MR 
(labeled E, I, K, 0, and R). 
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Climate change has a net negative economic impact on 76-88% of fanners with 
the future agricultural production system. The gains in the future under different 
GCMs range between 9.5 and 10% of mean farm returns (bottom left Fig. 11); 
losses range from 19 to 29% of mean farm returns). The losses under HADGCM- 
ES2 (GCM: K) are the highest. Net impact due to climate change would be negative 
for all GCM and losses shall range from 8-19% of net farm returns). 

Poverty rate in the future production system without climate change is 27.8%. 
However poverty rates would increase with climate change by 12-26%. Poverty 
rates and per capita incomes under Question 2 have decreased in comparison to 
Question 1 (without climate change). The increase in off-farm income by 30% as 
anticipated under the RAPs has caused this decrease. The decrease in household 
sizes (by 16%) has also contributed to bringing down poverty rates. 


Sensitivity analysis 

To address the sensitivity of the TOA-MD results to the choice of future trends 
in the RAPs, we have redone the analysis with a different set of growth factors 
obtained mostly from the global IMPACT economic model (Group 2 in Table 5). 
Between Group 1 and 2, the trends for global prices, household size, farm area, 
and non-agricultural income were identical. There were significant differences in 
the trends for yield, local prices, and cost of production with the values m Group 2 
being approximately about 15-20% greater. The results are shown in Table 7 and 
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Fig. 12. Variation of gains, losses, and net impacts as a percentage of net farm returns for Question 
2 utilizing Group 1 and 2 trends. 
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Fig. 1 2. Using the Group 2 trend values, the change in projected returns under climate ' : 
change were actually more positive for four GCMs in contrast to the results with > 
Group 1. Similarly the percentage change mpercapita income was favorable for f 0Ur 
GCMs under climate change, again in contrast to the results with trend factors from 
Group 1. These results show that the trend factors for yield, local prices and cost of ■ 
production are critical. They also show the high sensitivity of the results to the input ? 
trend factors. Moreover the percentage of vulnerable farms {% of losers) is much ■ 
ower under the Group 2 analysis compared to Group 1 . This provides a confirmation ' 
that the choice of trend factors from global economic models and RAPs are likely 
to be one of the vital key elements in the process of designing adaptation strategies J 

Conclusions and Discussion -"V 

A sustained program of research under AgMIP has led to the implementation of? 
clunate, crop, and economic models of rice agricultural systems in Sri Lanka. ' ri 

• >.; 

Conclusions • 

, 

Climate 

The climate projections show a clear increase in temperature in the mid-20th centuiy * 
for a high greenhouse gas concentration pathway (RCP8.5); these models show a 4 
tendency towards a wetter future but there is scatter in these results. The results were 
largely consistent across five locations in Sri Lanka and the details are provided for % ' 
data from the northwestern Province. The changes in the model output are described t 

Maha and Yala seasons "ext- 'Fte increase in average temperature for 4 
2040-2070 from that of 1980-2010 ranged from 1 . 1-2.4°C for Maha and 1.5-2.8°C 4 ‘ 
for Yala seasons. The rainfall shows a slight increase for a preponderance of the 20 '4 
GCMs in the CMIP5 archive, with the highest increase for October-November. £" 
During ; Maha, the average rainfall rises from a present average of 6 mm/day up to a 
future 9.5 mm/day for 19 of the models while it drops by 0. 1 mm/day for one model 4 
1 he magnitude of the rise was statistically significant for 16 of these models. During ? 
the Yala season, half of the models showed a decline in average rainfall of up to 1 ‘‘ ' 
mm/day from the normal of 2 mm/day, including four of five selected GCM outputs ; 
and 10 of the other models showed an increase — this increase was statistically 

significant for these ten models. 


Genetic coefficients for the DSSAT and APSIM rice models, which were generated 
with Set 1 data, are presented here. Genetic coefficients obtained independently with 
a subset of data are presented in Zubair et al. (2014) and are similar. 
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Simulations were undertaken for farmer fields with Set 1 and Set 2 of data with 
DSSAT and APSIM rice models and the DSS AT results are presented for both sets. 
These simulations showed an adequate match of the probability of exceedance in the 
case of DSSAT for both sets while the estimates from the APSIM models were lower. 

Under the projected climate change, the simulations using DSSAT with climate 
change scenarios projected relative yields that were lower for the Maha season (by 
6-19% for Set 1 and 6-15% for Set 2) and for the Yala season (24-41% for Set 1 
and 25-37% for Set 2). 

Crop simulations from DSSAT and APSIM with Set 1 showed contrasting results. 
The comparison of simulated and observed yields for the farmer fields in Set 1 
showed a lower RMSE with DSSAT than APSIM for Maha but not for Yala. Simu- 
lated yields showed yield drops for the future- 12-19%* for Maha and 27-41 % for 
Yala. 

The adverse effects were more pronounced for the Yala season. The GCM with 
the highest temperature increase — HadGEM2-ES — showed the highest yield 
drop — 41 % for Yala. 

Integrated assessment 

The TOA-MD economic model was used to assess the socio-economic impacts on 
paddy farmers due to climate change. The projected climate change brings losses for 
most farmers (ranging from 75-85% losers depending on the GCMs) for Migalewa. 
TOA-MD results for the GCM with the highest temperature increase (HadGEM2- 
ES) showed the lowest yields and highest economic impacts with 85% losers. 

The impact of climate change on the future production system was assessed 
(Question 2) assuming that future production system will differ from current 
production system due to continuation of current development trends in the agricul- 
ture sector. Considering such trends along with climate change leads to a simulated 
77-88% losers in Migalewa. 

The results were sensitive to the choices in parameter trends obtained from the 
RAPs and the ones obtained from IMPACT. This has important implications for 
the design and implementation of policies and adaptation strategies based on the 
analysis of possible future impacts of climate change. More work is needed to 
understand better the effects of technology and price trends as a result of global 
economic models or other projections. 


Discussion 

Climate 

The future projections for temperature reported here are consistent with past IPCC 
projections for South Asia (IPCC, 2007). The IPCC Fifth Assessment Report (AR5) 


'■ 
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anticipates an increase in temperature of 2.4 to 2.6°C in the global average for ' 
RCP2.6 for 2046-2065 and an increase in annual rainfall of up to 10% for 208 1 
2100 (IPCC, 2013). 

There has been a multiplicity of downscaled climate projections for Sri Lanka ' 
based on assessments prior to AR5 which was reviewed by Eriy agama and Smakthin 
(2010). These projections agree reasonably among themselves on projections for 
temperature and disagree widely among themselves for projections of precipitation. 

The contribution of this work is to undertake assessments with the entire ensemble of % 
CMEP5 models and present the ensemble of results, which also includes uncertainty 5 
among model implementations. 3 t % 

While this work provides for uncertainty among the climate models, we do need ¥' 
to address the larger issue of how well the models perform in this region, the small 
spatial scales that we consider, and the robustness of the models as they take on V' 
different regimes in the climate system. The climate models as implemented also do V. 
not provide a basis for future changes in extreme events. These are issues that could = ' 
be better answered in the future and should be considered in the interpretation of % '■ 
implications for policymaking. 

In the recent decades, there has been a tendency towards drying conditions in ii ' 
the study region and in Sri Lanka generally. Such a trend has not been captured 
by models in the past This phenomenon of near-term climate change needs to %■ ' 
be addressed in the future. Thus the paradox of the recent drying trend alongside 

projections of wettening needs to be addressed. V‘ . 


Crop 


We have assessed the sensitivity of rice yield with temperature, CO2, and r ain fo il by ’ 
varying each of these parameters individually in the anticipated range of changes 
with a DSSAT crop model implementation with Set 1 (Zubair et al„ 2014) The 
model shows very little sensitivity to changes in rainfall by, and mode’st sensitivity 
to C0 2 , but high sensitivity to changes in temperature. These results are consistent 
with past work that shows that rice yields are expected to increase with elevated CO- 
levels and decline with increases in temperature (Wassmann and Dobermann, 200?) 
but provide specific quantitative details for our model implementation in Zubair 
etal. (2014). 

The effects of pests and diseases on paddy production can override some of the 
projections and thus need to be taken account of in the future. The role of rising 
salinity will have to be considered in the coastal areas that are vulnerable to salt 
water intrusion into groundwater associated with sea level rise. 

There is a need to implement modules and test their implementation to take better 
account of water stress in the rice crop models as we use these models for relatively 
drier areas with limited irrigation in Sri Lanka. 
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There is a range of climates even within the rice-growing areas of Sri Lanka, 
which includes warmer, colder, drier, and wetter climates, and other soil types and 
associated agricultural practices for which the crop models should be further tested. 
One could also take advantage of practices in locations that may have a present 
climate similar to a projected future climate elsewhere. 

Integrated assessment 

In the RAPs that we used, die future agricultural practices help improve rice produc- 
tivity trends and help, mitigate the negative impacts of climate change on agricultural 
systems. Any RAP is only representative and other choices should be considered. 
The results from the TOA-MD model are highly dependent on the crop model sim- 
ulations and on model implementation. 

Any biases in the crop simulations would affect the TOA-MD results signif- 
icantly. The samples used to parameterize the model are fairly small under this 
study. In the future, the crop and economic models shall be implemented for other 
sites where we already have survey results. 

Our analysis points to the likelihood of inadequate income stability of farming 
and this can cause farmers to look for other employment options. Increase in off- 
farm employment under the RAP could have been one of the factors leading to this 
outcome. However this could pose a serious threat to rice production and will need 
to be addressed. 

An aggregate study (Seo et al., 2005) to assess climate change impacts on agri- 
culture in Sri Lanka show that impacts due to climate change could vary from -1 1 
billion rupees to +39 billion rupees nationally depending on the climate scenario 
(Seo et al, 2005). This work has provided for economic impacts on a much finer 
scale, taking better account of local conditions. Although one may say that there is 
general agreement with our work, the state of such studies are preliminary and highly 
dependent on future projections, methodology, and future agricultural pathways. 


Adaptation measures 

The adaptation strategies that were described in the first section of this chapter are 
aimed to deal with three important issues. Coping with extreme risks for farmers, 
varietal improvement, and coping with water scarcity. An extension to the current 
analysis will be to assess the impacts of those adaptation strategies on the agricultural 
production system of Migalawa. 

Farmer risk reduction and risk transfer 

Both our surveys of farmers and work with the economic analysis highlighted 
the need for better management of risks that fanners face. Better weather- and 
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climate-based hazard predictions, predictions for water availability, market-based 
risk predictions, the use of risk transfer mechanisms such as insurance and farming 
cooperatives, and better implementation of guaranteed price schemes shall have to- 
be prioritized if there is a policy priority for sustaining the number of rice farmers 


Changes in varietal use and varietal improvement •• 

There is significant variation of climate in Sri Lanka and agricultural practices. Thus ' 
as climate changes it may still be possible that varieties adapted to other regions 
become better suited for a given region. Such an analysis of climate shall point 
to some options for varietal adoption. Further analysis of data from varietal trials 
may help. In addition, there is an ongoing program of varietal improvement in Sn •" 
Lanka which may be critical as climate changes. The adoption of varieties that 
have enhanced flood, drought, salinity, and disease tolerance as the climate and soil : 
conditions change may be critical to tolerant rice varieties. 

s yy 


Water management ; :T; 


Our work here has been with assumptions of perfect irrigation. However as seen 
in Kadawaramulla, drought and irrigation failures are critical. Moreover, there is 
increasing likelihood of diversion of water for urban use and other purposes (Zubau 
et a/., 2014). Thus greater attention shall have to be paid to assessing the risk of water 
scarcity in the future in agricultural areas, the risks in ensuring adequate irrigation, 
and the use of water conservation techniques. 
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